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Abstract: We introduce a new family of multivariate wavelets which are ob- 
tained by "polyharmonic subdivision" . They generalize directly the original compactly 
supported Daubechies wavelets. 
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1 Introduction 

We consider new multivariate polyharmonic Daubechies type wavelets which 
are called " polyharmonic subdivision wavelets" . They have been recently introduced 
in the paper [5]. They are obtained by means of a procedure called "polyharmonic 
subdivision" which is a generalization of the classical one-dimensional subdivision 
scheme of Deslauriers-Dubuc [4] which is the original source for the first compactly 
supported wavelets of Daubechies in 1988, cf. [3 . This new family of polyharmonic 
wavelets is the second representative of the Polyharmonic Wavelet Analysis following 
the "polyspline wavelets" which have been introduced in the monograph [6]. 

An important feature of these newly-born wavelets is that they are a nice gener- 
alization of the one-dimensional wavelets of Daubechies: they form an orthonormal 
family, enjoy nice non-stationary "refinement operator" equations, and have compact 
filters. In addition to that they have elongated supports. Let us remind that a major 
drawback of the one-dimensional spline wavelets of Ch. Chui is that they do not have 
finite filters, and respectively, the polyspline wavelets of [6] do not have finite filters. 

2 Construction of fundamental function $m for exponential 
polynomials subdivision 

The whole construction of the Daubechies type wavelets passes via the construction 
of the so-called fundamental function of subdivision, cf. 1 . In the present case we 
will work with non-stationary subdivision and we have a family of such functions 
4>m for all m G Z which satisfy the refinement equations (two-scale relations) 
given by 




for all t e M. 



(1) 



We define the non-stationary subdivision symbol by putting 




(2) 



We are interested in special subdivision processes arising through the solutions of 
Ordinary Differential Equations. We assume that we are given a number of frequen- 
cies < Ai < A2 < ... < Ap and put for the frequency vector (with repetitions) 

A = {Ai, A2, Aat} U Ai, — A2, — Aat} . 

We consider the space of C°° solutions of the ODE 

j=i ^ / 

Let us recall a simple fact from ODEs: in the case of different Aj 's the space of all 
C°° solutions in Q is spanned by the set |e'^^* : j = 1,2, . In the case of s 
coinciding indices A.; = A^+i = ... ~ Xi+s-i we have that the solution set contains 
the functions {t^e^'* : £ = 0, 1, s} . 

Let us proceed to the construction of the subdivision symbols. We put 



k + l 



We define the following Laurent polynomial 

N 



d{z) :=#-l(z) ^^-Q (^ + ^.)(^"'+^.-) 



i=i 

and 



P,.,:.pl'IW:.n(l-^.) 



(4) 



They satisfy the equality 

d (e*") = P (sin^ I) for all w e R; (5) 

cf. [7]. We will often drop the dependence on the upper index in d, a, P and the 
other functions and symbols. 

An important step for construction of the subdivision coefhcients a^™' is the 
application of the Bezout theorem: 

Proposition 1 There exists a unique polynomial Q with real coefficients of degree 
N — I such that 

P (x) Q{x) + P{1- x)Q{l- x) = 1 



Q{x)>0 for xe {0,1). 
We define now the trigonometric polynomial b (z) — bl*^! (z) by putting 

b (e^-) ^ Q (sin^ I) . 

We finally define the symmetric Laurent polynomial a (z) by putting 

a(z) flW (z) := 2d(z)6(z) for z £ C \ {0} . (6) 

The following proposition is important for the application of the Riesz lemma to 
a (z) and construction of the Wavelet Analysis, cf. [7], [5]. 
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Proposition 2 The polynomial a (z) defined in (0) satisfies 

2N-1 

a{z)= ajz^ 

j=-2N + l 

with aj — Q-j — aj and 

a{z) > for all \z\ — 1. 

The following fundamental result shows that the symbols a (z) are the non- 
stationary subdivision symbols for symmetric set of frequencies A, cf. [5]. 

Theorem 3 For every exponential polynomial, i. e. for every solution to the equation 

Lfit)--=Ili^~^^)fW=0 (7) 

we put 

Then f is reproduced by means of interpolatory subdivision^ i.e. 

CO 

fj^'^ E 4-2, f J for all fez (8) 

j=~oo 

= f^ for all J e Z, 

For every m G Z the fundamental function of subdivision <I>„i (t) is a continuous 
function obtained throught the subdivision process (0), where one starts from /j* = 5j 
for j £ Z (here 6j is the Kronecker symbol), i.e. we put $m (2^) = <5j, and 
satisfies the refinement equation (QJ). 

Having in hand the functions $,„ and their refinement symbols a[™l we may follow 
the usual scheme for construction of father and mother wavelets which has been used 
by Daubcchies, cf. [3], [1]. The following fundamental result has been proved in [5]. 

Theorem 4 There exists a polynomial g (z) ~ X^jez that it is the "square 

root" of 2a (z) , i.e. 

a{e^') = l\g{e^')\' (9) 

For every m £ Z there exists a compactly supported function (pm (t) which satisfies 
the refinement equation 

{t) = gjfm+l (2t - j) , (10) 

j 

and the family {fm {t — orthonormal. (These are the non- stationary father 

wavelets.) The functions 

i^m it) = J2 (-1)' 5l-J</'m+l {2t - j) (11) 

are the mother wavelets; the family {ijjm {t — j)}jgz is orthonormal and the family 
{ipm [t — j)}m j£Z forms an orthonormal basis of L2 (M) . 
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2.1 The polyharmonic case 

For the polyharmonic subdivision we will work with very special ODEs defined by 
:= {(P /dt^ ~ 'f^) which are the Fourier transform of the polyharmonic operator 
A^. For a fixed constant ^ > we put 

A:=(-^,-e,...,-e,^,C,-,OeM'^ (12) 
i.e. \j = ^, for j = 1, 2, A^. Now for fixed £,>Q and fc G Z we define the polynomial 

d (z) := (z) := rfW (z) := "^"^ """^ for z S C; (13) 

(1 + a;o) 

here we put xq := e~^l'^ ^ . For the sake of simplicity we will very often drop the 
dependence on A; and ^. By ([S]) we have d (e*'^) — P (sin^ ^) where 

P(:r)= l--i^J =(1-^?^^, (14) 
V (1 + xo) ) 

and we have put 

Jfc],5 4x0 2 



77 ry^ 



(1 + xo)' l + cosh(e/2fe+i)- 

Then following Proposition [T] we have to find the polynomial solution Q to the equa- 
tion 

P (x) Q (x) + g (1 - x) P (1 - x) = 1 
where Q has degree < — 1. 

Remark 5 Let us recall that the polynomial Q in the classical case, cf. e.g. J3]/, p. 
195, satisfies condition 

{l-yfQ{y) + y''Q{l-y) = l. 

The lowest degree solution polynomial Q will be called Daubechies ' polynomial and 
we put 

i?^(x):=Er^'"')^^- (15) 
j=o V / 

(Note that in J^I and ^ the notation used is Ppj !) 

It is amazing that it is possible to solve the problem in Proposition [T] explicitly. 

Proposition 6 Let A = (— ^, — ^, — .-.,£,) 6 K^^. Then for the correspond- 
ing polynomial P (x) = (1 — rjx) , the polynomial Q of degree N — 1 defined by 

fN + j-l\ (l-77(l-x))^ 

j=o 

solves the equation 

F(x)g(x) + p(i-x)g(i-x) = i. (17) 

Hence, 

g(x) = (2-,)-^P^(i^|il-^). (18) 



g (x) = g^« (X) = (2 - ,)-- ^ - - ^ - M (16) 
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Hence, we find the trigonometric polynomial fct'^l (z) by putting 

foW (z) := 6W'« (e^-) gW^« (sin^ |) (19) 



where we recall the notations 

o a; 1 — cos uj 1 



X = sm — 



2 2 2 4' 

Finally, we obtain the subdivision symbol a^'^l (z) by putting 

(z) a['''l^« (z) := 2dW^« (z) 6^'^ (z) . (20) 

Now by Theorem|4]we find the "square root" of the symbol a^'^l (z) . This means that 
we have to take separately the "square root" of the Laurent polynomials d}''^ (z) and 
6^ (z) . The " square root" of d^'^l (z) is obvious; taking the "square root" of b^'^^ (z) 
needs taking the "square root" of the polynomial Q. 

3 Algorithm for finding the square root of the polynomials Q 

For the algorithmic aspects of taking the "square root" of the polynomial Q it will 
be important to describe the polynomial Q through the zeros of the Daubechies' 
polynomial Rn in (llSp . 



Proposition 7 Let the zeros of the Daubechies' polynomial U5]) be Cj', i. 

'^fN + j-l\ , (27V -2)! rj. 

Then the polynomial Q as determined by I116\) is given by 

N~l 



e. 



cf (2 - 77) + 77 - 1 



((^-1)!) ,=1 

where 



By formula (113^ we have the representation 



(z) 



(z + xo) 



for z 



hence, we take the trigonometric polynomial 

Mr (z) := ii±^ (21) 
(1 + 2:0) 

as its "square root", i.e. d^^^ (z) = |Mi (z)|^ for |z| = 1. Further, we have to take care 
of the "square root" of the polynomial 6^'^' (z) . Thus we have to find the polynomial 
M2 of degree < iV — 1 such that 



i0(W|). (22) 



which may be obtained by using the roots of the Daubechies polynomials. 
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Remark 8 Let the polynomial Q have the zeros Cj as in Proposition and let us 
put 

We see that Q (sin^ ^) ~ Q (cosw) for some polynomial Q and cj are the zeros of Q. 
Hence, we may apply the algorithm for the Riesz representation of Q, see e.g. JTI, p. 
197- 198. 

Thus we obtain finally for every integer m > and ^ € Z" the representation 

a^^U^\{z)^\\M,{z)Nh{z)\\ (23) 

and the family of functions 

M [z) := M[™1 (z) -.^ M[™1'« (z) := Mi (z) A/a (z) (24) 

represents the refinement masks for the family of scaling functions (father wavelets) 
{ipm (0}m>o which the functions $„i are autocorrelation functions. 

Remark 9 Note that the above factorization has been found in the special case ^ = 
by Daubechies in jSj, p. 266; the coefficients of the "square root" polynomial for 
N — 2. .10 are in table 6.1 in J^. A detailed discussion of more efficient meth- 
ods for choosing the proper polynomial M2 (z) is available in Strang-Nguyen p. 
157, in chapter 5.4 on Spectral factorization. The factorization of the Daubechies' 
polynomial Rn (y) is discussed in Burrus J^, on p. 78 and the Matlab program is 
[hn,hinJ—daub(N) in Appendix C. They work with the zeros of the polynomial Rn 
and provide a number of manipulations for finding a more stable factorization. 
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